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Simplidal and Cellular Trees 


Art M. Duval, Caroline J. Klivans, and Jeremy L. Martin 


1 Introduction 

How can redundancy be eliminated from a network, and what can be said about the 
resulting substructures? In graph-theoretic terms, these substructures are spanning 
trees. The collection of all spanning trees of a graph has the structure of a matroid ba¬ 
sis system; this observation connects trees to algebraic combinatorics and explains 
why many graph algorithms can be made computationally efficient. The number of 
spanning trees of a graph measures its complexity as a network, and there are clas¬ 
sical and efficient linear-algebraic tools for calculating this number, as well as for 
enumerating trees more finely. The algebra of trees is key in studying dynamical 
systems on graphs, notably the abelian sandpile model (known in other forms as the 
chip-firing game or dollar game), whose possible states are encoded by a group of 
size equal to the complexity of the underlying graph. 

This chapter is about the more recent theory of spanning trees in cell complexes, 
which are natural higher-dimensional analogues of graphs. The theory for cell com¬ 
plexes parallels the graph-theoretical version in many ways, including the connec¬ 
tion to matroid theory. However, higher-dimensional spaces can have much richer 
topology, which complicates the algebraic and enumerative parts of the story in a 
very concrete way. It turns out that the simple number of spanning trees of a cell 
complex is not a well-behaved invariant, and it is better to account for topological 
complexity in the enumeration. On the other hand, many essential algebraic tools 
for working with spanning trees of a graph do extend well to arbitrary dimension. 
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Sections 1 and 2 give an overview of the subject, including the history of the 
subject of higher-dimensional trees and emphasizing the geometry and topology of 
trees. As far as possible, we have attempted to make these sections self-contained 
and accessible to readers not familiar with the machinery of algebraic topology, al¬ 
though a certain amount of technical material is unavoidable. The next two sections 
are more technical: Section 3 describes the various extensions of the matrix-tree 
theorem from graphs to cell complexes, and Section 4 examines individual classes 
of cell complexes whose tree enumerators are understood. Finally, Section 5 covers 
new directions, including extensions of algebraic graph theory and connections to 
critical groups, cuts and flows, simplicial decompositions, and matroids. 


1.1 Higher-dimensional trees 

The higher-dimensional trees considered here were first studied by Bolker [Bol76], 
in the context of transportation polytopes, which are solution sets to certain combi¬ 
natorial optimization problems. The vertices of transportation polytopes are span¬ 
ning trees of complete colorful complexes, which are higher-dimensional general¬ 
izations of complete bipartite graphs. Bolker pinpointed the basic homological con¬ 
ditions on spanning trees (essentially Definition 2.6 below). Kalai [Kal83] consid¬ 
ered spanning trees of complete simplicial complexes, which generalize complete 
graphs, and obtained a beautiful enumerative formula generalizing Cayley’s for¬ 
mula. Soon after, Adin [Adi92] gave a corresponding formula for counting span¬ 
ning trees in complete colorful complexes. More recently, many authors, including 
Duval, Klivans and Martin [DKM09, DKMll, DKM15], Lyons [Lyo09] and Pe- 
tersson [Pet09], independently observed that Bolker and Kalai’s definitions could be 
extended to general cell complexes to develop a broad theory of higher-dimensional 
trees. 

Algebraically speaking, one can define a spanning tree of a suitably connected 
cell complex to be any subcomplex generated by faces corresponding to a linearly 
independent set of columns of a certain matrix (see Definition 2.5). However, this 
definition does not really reflect the topology of higher-dimensional trees. By anal¬ 
ogy, a spanning tree T of a connected graph G can be defined as a collection of 
edges which correspond to a column basis of the oriented incidence matrix of G. 
But this definition does not directly capture what a tree looks like, and the most 
familiar definitions are more topological and combinatorial in nature: 

• r is connected and acyclic; 

• r is connected and jTl =n—l (where n is the number of vertices); 

• r is acyclic and \ T\ = n — 1-, 

• r is a maximal acyclic subgraph of G; 

• r is a minimal connected subgraph of G. 
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Similarly, for a c/-dimensional cell complex X with codimension-1 Betti number 
Pd-i{X) = 0 (the analogue of connectedness for graphs), a subcomplex T C X is a 
spanning tree if it satisfies any of the following equivalent properties: 

• j5diT)=j5d-i{T)=0-, 

. l5d{T)=0and\Td\ = \Xd\-Mxy^ 

• j5d-i{T)=Oand\Td\ = \Xd\-l5dixy, 

• r is maximal among the spanning subcomplexes of X with PdiT) =0; 

• r is minimal among the spanning subcomplexes of X with l5d-i{T) = 0. 

A more complete list of characterizations, including the more general notion of 
spanning forests when l3d-i{X) > 0, is given in Proposition 2.13. Definitions and 
notation for the terms above are also given in Section 2. 

Many cellular spanning trees “look like” trees. For instance, a 2-dimensional 
simplicial complex whose dual graph forms a tree automatically satisfies the con¬ 
ditions above; see the first complex in Figure 1. On the other hand, many higher¬ 
dimensional trees look less familiar. For instance, if X is any cellular sphere (such 
as the boundary of a polytope), then its spanning trees are the contractible subcom¬ 
plexes formed by removing a single maximal face; this is the higher-dimensional 
generalization of the fact that the spanning trees of a cycle graph are the paths 
formed from it by deleting one edge. Of particular interest (and the source of much 
difficulty) is the existence of trees with torsion homology, which can occur in any 
dimension > 2. Already in dimension 2, trees can have finite non-trivial homology. 
The smallest simplicial example is the standard triangulation of the real projective 
plane which arises as a spanning tree of the 2-dimensional skeleton of the 
6-vertex simplex. See Examples 2.11 and 2.12. 


1.2 Tree numbers 

The subject of enumerating spanning trees traditionally dates back to Kirchhoff’s 
1847 work on electrical circuits [Kir47]. Excellent general references with more 
historical details include Moon [Moo70] and Stanley [Sta99, pp. 54-69]. In modern 
terms, let G be a connected graph on n vertices whose Laplacian matrix L = L{G) 
has nonzero eigenvalues ,..., i, and let ^(G) denote its set of spanning trees. 
Then the matrix-tree theorem (or Kirchhoff’s theorem) states that 

t(G) = detZ = - hxiL ( 1 ) 

n 

where t(G) = and L denotes the reduced Laplacian obtained from L by 

deleting the row and column (for any i). Thus tree enumeration is closely tied in 
with the broad subject of spectral graph theory; see, e.g., [CDS80, Chu97]. A most 
important special case is the formula 


n—2 


x{Kn) = ri 


( 2 ) 
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for the number of spanning trees of the complete graph Kn, known as Cayley’s for¬ 
mula, although it was in fact first observed by Sylvester [Syl57] and first proved, in 
an equivalent form, by Borchardt [Bor60]. Cayley’s formula can be deduced from 
the matrix-tree theorem by finding an explicit eigenbasis for L{K„). Independent 
combinatorial proofs include the well-known Priifer code [PriilS] and Joyal’s bijec- 
tion using permutations [JoySl, pp. 15-16]; see also [Sta99, pp. 25-28]. 

There are other families of graphs with nice spanning tree counts which can be 
obtained either algebraically (using the matrix-tree theorem) or combinatorially (via 
a bijection). For example, the number of spanning trees of the complete bipartite 
graph is 

(3) 

as discovered by several authors [Aus60, SH59, FS58, Sco62], and arises in the 
transportation problem of combinatorial optimization [KW68]. 

Moving to higher dimensions, in Bolker’s original work, he observed that it is 
“tempting on numerological grounds” to generalize the well-known formulas (2) 
and (3), but also was aware that the obvious generalizations were incorrect because 
of torsion trees, i.e., c/-dimensional trees T with ;Z) finite but nontrivial. 

Kalai [Kal83] was the first to pinpoint the role played by torsion. He showed that 
the correct analogue of the tree count for a cZ-dimensional simplicial complex X is 
not simply the cardinality of the set fl'iX) of simplicial spanning trees, but rather 
the quantity 

T(X)= Y. \Hd-i{T)\^ (4) 

TeS{x) 

where denotes the {d — 1)-dimensional reduced homology group of T 

over Z. The summands in (4) arise naturally in the expansion of the combi¬ 
natorial Laplacian matrix of X using the Binet-Cauchy formula. Note that tor¬ 
sion has no effect on spanning tree enumeration in graphs; every topological 
space has torsion-free 0* homology group, so when 2f is a graph the right-hand 
side of (4) reduces to the number of spanning trees. (On the other hand, torsion 
does arise in enumerative problems equivalent or related to tree enumeration; see, 
e.g., [MR05, BH06, BBGM14, DM12].) Additional evidence that T(2f) is the right 
way to count trees is Kalai’s striking generalization of Cayley’s formula: 

(5) 

where ^ denotes the cZ-dimensional skeleton of an n-vertex simplex. Subse¬ 
quently, Adin [Adi92] used the same ideas to give a formula for the torsion-weighted 
count of a complete colorful complex, generalizing (3). Both formulas were pre¬ 
cisely those predicted by Bolker. Kalai also gave a weighted version of (5) enumer¬ 
ating spanning trees by their degree sequences [Kal83, Thm. 3'] and showed that 
torsion accounts for most of the quantity T(2f): as Kalai put it, spanning trees of 
K„j are “on the average, far from being [Z-]acyclic” [Kal83, §5]. 

The torsion-weighted spanning tree count of a cell complex X‘^ can be computed 
from its combinatorial Laplacian operator L. This essential fact, the cellular matrix- 
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tree theorem, has appeared in various forms in [DKM09, DKMll, Lyo09, Pet09, 
CCK12], all building on the ideas originating in [Kal83] and [Adi92]. We will de¬ 
scribe the different guises of this theorem, and their interconnections, in Section 3 
below. As in the graphical case, t{X) may be expressed via a determinant or via 
eigenvalues. Specifically, for a cZ-dimensional complex X, the torsion-weighted tree 
count of c/-dimensional trees is given by 


Xd{X) = c-detZ = c' 


Td-liX) 


where c and c' are certain correction factors arising from torsion homology; L is the 
reduced Laplacian formed from L by removing the rows and columns corresponding 
to a codimension-1 spanning tree; and the A, are the non-zero eigenvalues of L. 

In practice, for many families of complexes of combinatorial interest (e.g., 
Cohen-Macaulay complexes), the homological correction factors c and c' disappear. 
In this case, the weighted tree count is given precisely by a determinant, as in Kalai’s 
enumeration of skeletons of simplices. In general, the eigenvalue formulation of the 
matrix-tree theorem is recursive; the count for c/-dimensional trees is expressed in 
terms of the count for {d — 1)-dimensional trees. In the case that c' = 1, the recursion 
yields the alternating-product expression 


Td{X) 


Ai■A3■As■ ■ ■ 
Aq•A2•A4• • • 


where A, is the product of the non-zero eigenvalues of the /th Laplacian. 

Laplacians of simplicial complexes first appeared in the work of Eckmann [Eck45] 
and Dodziuk and Patodi [DP76]. They were studied as combinatorial objects by 
Eriedman [Eri98] and Eriedman and Hanlon [EH98], who in particular considered 
the phenomenon of integral Laplacian eigenvalues. Subsequently, many combinato- 
rially significant classes of complexes have been shown to have all integer Laplacian 
eigenvalues, such as chessboard complexes [EH98], matroid complexes [KRSOO], 
shifted complexes [DR02], and hypercubes [DKMll]. The strong connection be¬ 
tween Laplacian eigenvalues and tree enumeration suggests that these families 
might support nice spanning tree counts with combinatorial interpretations. Indeed, 
there are many successes in this direction including enumerators for shifted com¬ 
plexes in terms of critical pairs [DKM09] and for matroid complexes in terms of 
flats and the j3-invariant [KL15a]. We will look at many of these families individu¬ 
ally in Section 4. 


1.3 Beyond basic definitions and enumeration 

The theory of trees has recently emerged as a natural aspect of the combinatorics of 
cellular spaces. We describe several new directions in Section 5. Many fundamental 
objects of algebraic graph theory have been generalized to higher dimension, in- 
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eluding critical groups and generalized sandpile/chip-firing models (§5.1) and cuts 
and flows (§5.2). Simplicial trees provide an inroad to studying theorems and con¬ 
jectures about decompositions of simplicial complexes (§5.3). Finally, the study of 
cellular matroids has led to new topological invariants on cellular manifolds and 
applications to probability (§5.4). 


2 Trees and forests: from graphs to cell complexes 

The definitions of higher-dimensional trees and forests may be approached from 
various directions, including linear algebra and matroid theory, but ultimately some 
topological machinery is necessary to appreciate the ideas fully. We will attempt to 
give an informal, self-contained description of those concepts that we need such as 
cell complexes and homology, that will be accessible to a reader willing to endure a 
little black-box algebra. For the precise details, the reader should consult a standard 
reference such as Hatcher [Hat02]. 


2.1 Trees and forests in graphs 

Let G = {y,E) be a finite graph, with |y| = n and iFj = m. For each edge e on 
vertices v and w, fix an orientation e = v —> w. The signed incidence matrix Bq is 
then defined as the n x m matrix with entries 

{ 1 if e = M —)■ V for some ufv, 

— 1 if e = V —)■ w for some w fv, 

0 otherwise. 

In particular, the column of d is zero if and only if e is a loop. We now recall 
several well-known equivalent definitions of trees and forests from graph theory. 

Definition 2.1. The graph G is a tree if it satisfies any of the following equivalent 
properties: 

(a) G is connected and acyclic. 

(b) G is connected and m = n—1. 

(c) G is acyclic and m = n—1. 

(d) Every pair of vertices in G is connected by exactly one path. 

(e) G is a maximal acyclic graph. 

(f) G is a minimal connected graph. 

(g) For any (hence every) orientation of G, the columns of dg are a basis for the 
space Mq = {V G M" | vi H-h v„ = 0}. 

Meanwhile, G is a forest if it satisfies any of the following equivalent properties: 
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(h) G is acyclic. 

(i) Every pair of vertices in G is connected by at most one path. 

(j) The columns of the incidence matrix Bq of G are linearly independent. 

(k) Every connected component of G is a tree. 

In order to extend the ideas of trees and forests, we first have to describe the 
objects that will serve as higher-dimensional analogues of graphs. 


2.2 Cell complexes 

A cell complex X is a topological space that is the disjoint union of cells, each 
homeomorphic to the interior of a k-dimensional disk for some k, attached together 
in a locally reasonable manner. In particular, the closure of any cell is a union of 
cells. The set of k-cells is denoted Xj^, and the k-skeleton is the union of all cells 
of dimension < k. A cell not contained in the closure of any other cell is called a 
facet, and a complex whose facets all have the same dimension is called pure. Every 
subset S of cells generates a subcomplex Xs = Uaes (We often blur the distinction 
between a set of cells and the subcomplex it generates.) The dimension of X is the 
largest dimension of one of its cells; the notation X‘^ indicates a cell complex of 
dimension d. In particular, a graph is exactly a cell complex of dimension 1. 

The combinatorial data about how cells are attached is recorded by a sequence 
of linear maps dj^ = di^{X) \ —)■ for Q <k < d. The map d\ is just the 

(signed) incidence matrix of the 1-skeleton graph, and the other maps are higher¬ 
dimensional analogues of the signed incidence matrix of a graph. Eor k>Q, each dk 
may be represented as a matrix with rows indexed by {k— l)-cells p and columns 
indexed by k-cells a; the entry dp a is an integer that records, roughly speaking, the 
number of times that a is wrapped around p, counted with orientation. If p is not 
contained in the closure of a, then dp a — 0. In particular, if some {k — I)-celI p is 
not in the closure of any k-cell, then the corresponding row of d is zero. We formally 
define (9o : —> M to be the map that sends every point to 1 e M. 

It is a standard fact that im(9i.+i C keri?^. for all k, and the quotient Hi^{X-,'R) = 
ker(9^/im(9i.+i is a topological invariant of X called its k* (reduced) homology 
group. The k* (reduced) Betti number is pi^lX) = dimK//i;(X;IR). (All homology 
groups and Betti numbers appearing here are reduced and we will use the notation 
Hk and in place of the usual and for ease of readability.) The 0* reduced 
Betti number is one less than the number of connected components, so in particular 
Po(X) = 0 if and only if X is connected. If G = (ViE) is a graph with c connected 
components, then j3i (G) = jE | — (jV | — c) is the size of the complement of a maximal 
spanning forest (sometimes called the cyclomatic number). 

Simplicial complexes are the cell complexes most familiar to combinatorialists. 
An (abstract) simplicial complex is a set A of subsets of [n] = {1,2,..., n}, such that 
whenever a € A and p Qo, then p G A. Thus a simplicial complex is completely 
determined by its set of facets. An abstract simplicial complex can be regarded as a 
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topological space by associating O’ = {t'l,..., /jt} € 4 'with the interior of the convex 
hull of ,... in K”, 'where ej is the standard basis vector. For readers 
familiar 'with hypergraphs, a simplicial complex is just a hypergraph in 'which every 
subset of a hyperedge is also a hyperedge. 

Simplicial complexes are regular cell complexes, 'which means that every attach¬ 
ing map is locally a homeomorphism; i.e., one-to-one. In particular, for a regular cell 
complex, every entry of every matrix dk is in {0,±1}. 

Example 2.2. A graph G is a 1-dimensional cell complex 'whose attaching map d\ 
is just the signed vertex-edge incidence matrix. As a cell complex, G is regular if G 
has no loops, and simplicial if in addition G has no parallel edges. 

Example 2.3. Gluing t'wo hollo'w tetrahedra along a common triangle 123 produces a 
cell complex, the equatorial bipyramid, 'with seven 2-cells (triangles), nine 1-cells 
(edges), and five 0-cells (vertices). This is an example of a simplicial complex, 'with 
facets {123,124,125,124,134,234,235} 'where 1,2,3 are the “equatorial” vertices 
and 4,5 are the “poles”. The collection of facets 'would often be considered a 3- 
uniform hypergraph; i.e., a hypergraph 'where all edges have size 3. 

Example 2.4. The real projective plane M.P'^ can be realized as a cell complex 'with 
one cell at of each dimension k for k G {0,1,2}. The 1-skeleton CJi U (Jo forms 
the line at infinity; the attaching map di is the matrix [0] because CJi is a loop. 
Mean'while, the attaching map 82 is the matrix [2], representing the fact that 02 
is 'wrapped t'wice around CJi. The smallest triangulation of M.P^ (i.e., the smallest 
simplicial complex homeomorphic to M.P^) can be obtained by identifying antipodal 
faces on an icosahedron. See Figure 3. 


2.3 Trees and forests in cell complexes 

We can no'w define the cellular analogues of trees and forests. 

Definition 2.5. A cZ-dimensional cell complex A is a cZ-forest if the columns of 
dd{X) are linearly independent. A spanning (/-forest of A is a (/-forest PCX such 
thatF<d-i =X<d-i. 

Here the term “spanning” refers to the requirement that the {d — 1)-skeleton of F 
equals that of A, just as a spanning subgraph of a graph G is a subgraph that includes 
all vertices. From a homological perspective, an equivalent condition for A to be a 
(/-forest is that //^(A;IR) = 0, because //(/(A<^;IR) = kevdj{X). 

Definition 2.6. A (/-dimensional cell complex A is a (/-tree if the columns of d^{X) 
are linearly independent and r?mkdd{X) — nullity 1 (A). A spanning (/-tree of A 
is a (/-tree PCX such that = A<^_i. 
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Equivalently, X is a tree if (X; M) = 1 (X; M) =0. 

Note that every cell complex has at least one spanning forest (since its bound¬ 
ary matrix has at least one column basis), but a spanning forest of X is a span¬ 
ning tree if and only if = 0. In general, the conditions = 0 

and = 0 should be regarded as the analogues of acyclicity and con¬ 

nectedness for a -dimensional cell complex. Thus Definitions 2.5 and 2.6 say that 
(f-forests and tZ-trees are respectively “c/-acyclic” and “t/-acyclic and cZ-connected” 
cell complexes. 

What do cellular trees and forests look like? A 1-dimensional forest is a graph 
that contains no cycles (subgraphs homeomorphic to a circle) and a 1-dimensional 
tree contains no cycles and is connected (i.e., every pair of vertices can be seen 
as the endpoints of some 1-dimensional path). Correspondingly, a 2-dimensional 
cell complex A is a forest if X contains no subcomplex homeomorphic to a compact 
oriented 2-manifold. Imagine a 2-dimensional cell complex as a network of bubbles: 
in order to make it into a forest, one must remove enough 2-cells to pop all the 
bubbles, without removing any lower-dimensional cells. Furthermore, A is a 2-tree 
if every 1-cycle is the boundary of some disk consisting of 2-cells. Definitions 2.5 
and 2.6 make sense even when d = 0: a 0-dimensional forest is any set of vertices, 
and a 0-dimensional tree is just a single vertex, which is a 0-spanning tree of A if 
and only if A<i is a connected graph. 

Example 2.7. Figure 1 shows three 2-dimensional simplicial trees. Complex (a) is a 
triangulation of the disk whose dual graph is a tree. Complex (b) is also a triangu¬ 
lation of the disk, but the dual graph forms a cycle, so it is a 2-tree with no “leaf”. 
Complex (c) is a contractible simplicial complex (hence a tree) with facets 123, 124, 
134, 125, 135. It is a spanning tree of the equatorial bipyramid. 





Fig. 1 Some two-dimensional trees. 


Example 2.8. If X‘^ is homeomorphic to a compact oriented manifold (such as a d- 
sphere) then PdiX) = 1, and the kernel of is one-dimensional, generated by a 
linear relation that involves all its columns nontrivially. Informally, A has exactly 
one bubble, and it can be popped by removing any facet. Deleting subsequent facets 
pokes extra holes in the surface. Hence any spanning subcomplex obtained from A 
by removing one facet (resp., one or more facets) is a maximal spanning forest 
(resp., a spanning forest). 
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Example 2.9. Consider the equatorial bipyramid B of Example 2.3. Topologically, it 
is a “double bubble”, recorded homologically by its second Betti number ^ 2 {B) = 2. 
In order to pop both bubbles, we can remove either (a) two triangles, one containing 
vertex 4 and one containing vertex 5, or (b) the middle triangle 123 plus any other 
triangle. Thus B has a total of (3 x 3) + 6 = 15 spanning trees, all of which are 
contractible simplicial complexes with 5 facets. See Figure 1(c) and Example 2.7 
for one example. 

Example 2.10. Figure 2 shows two 2-dimensional simplicial forests, which triangu¬ 
late the annulus and the Mobius band. In both cases p 2 = 0 but j3i = 1 > 0 (since 
neither complex is simply connected), so each complex has itself as a spanning 2- 
forest, but has no spanning tree — just like a non-connected acyclic graph. On the 
other hand, unlike a 1-dimensional forest, neither complex is the disjoint union of 
2-trees. 



The connectedness condition (2f ;IR) = 0 does not preclude the existence of 
torsion. Suppose that we regard the boundary matrix as a map on Z-modules 
. jXi, (rather than on M-vector spaces), and accordingly work with the 

integer homology groups 74(2f;Z) = df /While working over Z does 

not change the ranks of the matrices, the condition rank(9(;(2f) = nullity i (2f) 
no longer implies that //(;_i(2f;Z) is the trivial group, merely that it infinite. This 
subtlety is concealed in the graph case d =\, since it is a fact that Ho{X-,Z) is free 
abelian for any space X, but is of central importance in higher dimension. 

Example 2.11. While every tree graph is a contractible topological space, this is 
not always true in higher dimension. The simplest example is the real projective 
plane. In the standard cell structure on M.P'^ described in Example 2.4, \md 2 = 2IR = 
ker(9i(2f). So Hi{RP^-,R) = 0 and RP^ is a 2-tree (hence the only spanning 2-tree 
of itself). On the other hand, Hi(RP^;Z) = Z/2Z ^ 0, so RP^ is not contractible. 
In the triangulation of shown in Figure 3, the cycle consisting of edges 12, 13, 
and 23 is a representative of the non-trivial torsion homology class. 

Example 2.12. Let K^ 2 be the 2-dimensional skeleton of the 6-vertex simplex, i.e., 
the 2-dimensional simplicial complex on vertex set V = [6] whose facets are the 
twenty 3-element subsets of V. Let T\ be the spanning subcomplex of K(, 2 whose 
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facets are the ten 3-element sets containing vi. Then T\ is contractible (it is a cone), 
hence is a simplicial spanning tree of K(, 2 with //i(7’i;Z) = 0. On the other hand, 
M.P^ can be triangulated with six vertices (as in Figure 3), and any such triangulation 
T 2 is also a spanning tree of K(, 2 , with Hi (72; Z) = Z/2Z. 



2 


3 


Fig. 3 The triangulation of the real projective plane is a 2-tree with non-trivial torsion. 


Just as deleting an edge from a graph either decrements the cyclomatic number 
or increments the number of components, deleting facets of a simplicial complex X‘^ 
either decrements or increments jSd- 1 . Thus a (7-forest is a subcomplex T cX ob¬ 
tained by removing enough facets so that j5d{T) = 0, and a maximal (7-forest is one 
in which no facets have been removed unnecessarily, so that j5d-i{T) = l5d-i{X). 
If in addition (Jf) = 0, then T is a spanning tree. Thus the vanishing of j5d{T) 
and j5d-i{T) are the higher-dimensional analogues of the graph-theoretic concepts 
of acyclicity and connectedness respectively. 

The number of facets of every maximal (7-forest TCXis \Td\ = \Xd\-MX), 
which is the analogue of the statement that a maximal spanning forest of a graph 
with n vertices and c components has n — c edges. 

This discussion leads to the following characterizations of cellular spanning 
forests and trees, which generalize most of the conditions of Definition 2.1. 

Proposition 2.13. Let X be a cell complex of dimension d and let T contain all 
faces of dimension {d — 1) or less ofX. Then T is a maximal spanning forest ofX if 
it satisfies any of the following equivalent properties: 

(a) pd{T)^0andfid-i(T)=l5d-iiX). 

(b) l5d-i(T) = l5d-i(X) and \Td\ = \Xd\-MX)- 

(c) lidiT)^0and\Td\ = \Xd\-MX)- 

(d) Every element ofk&xdd-i{X,T) is the boundary of exactly one d-chain in T. 

(e) T is maximal among the spanning subcomplexes ofX with ^d{T) = 0. 

(f) T is minimal among the spanning subcomplexes ofX with pd-iiT) = fid-iiX). 

(g) The columns of dd{T) are a vector space basis for colspace(9^(Jf). 

Meanwhile, T is a spanning forest (not necessarily maximal) if 

(h) MT)=0. 

(i) Every element ofk&xdd-i{X,T) is the boundary of at most one d-chain in T. 
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(j) The columns of d^{T) are linearly independent. 

A maximal spanning forest is a spanning tree when = 0. In this case 

ksrdd-i{X,T) becomes ksrdd-\{X) in (d) and (i). 

When d=\, these properties reduce to the familiar equivalent definitions of trees 
and forests given in Definition 2.1, letter by letter. One condition that does not carry 
over to higher dimension is (k) of Definition 2.1, which defines a forest as a graph 
in which every connected component is a tree, it relies on the fact that the incidence 
matrix d\{G) of a graph G breaks up as a block sum of the incidence matrices 
of its components, of which there are j3o(G) + 1. On the other hand, there is no 
canonical way to decompose a t/-dimensional cell complex X, or the matrix dd{X), 
into j5d-i{X) + 1 pieces. Hence, c/-forests are not simply disjoint unions of (f-trees, 
as in the annulus and Mobius band of Example 2.10. 

Remark 2.14. If E C is a spanning tree, then there is a surjection of finite groups 
Hd-\{T','Z,) —> In particular, if X has a Z-acyclic spanning tree then 

Hd-i{X',1^) = 0, but the converse does not hold. For example, consider a cell com¬ 
plex with top boundary map [p q\ for any relatively prime integers p,q> 1. 


2.4 Rooted trees and forests 

What is the higher-dimensional analogue of a rooted tree or forest? In dimension 1, 
a root of a tree graph is a vertex, while a root of a forest graph is a choice of a vertex 
from each connected component. As before, in order to generalize these concepts to 
arbitrary dimension, we first state them in terms of linear algebra. It is a standard fact 
that the null vectors of a graph G are just the characteristic vectors of its connected 
components; therefore, a root of a maximal spanning forest of G can be regarded as 
the complement of a column basis of Bq, which motivates the following definition. 

Definition 2.15. Let X be a c/-dimensional cell complex. A rooted spanning forest 
of is a pair {F,R), where E is a spanning forest of 2f; /? is a (c/ — 1 )-dimensional 
subcomplex of X such that R<d -2 = X<d- 2 \ dx\R.F is square and nonsingular, 
where is the submatrix of dd{X) with rows X\R and columns F. 

Figure 4 shows rootings of some of the 2-trees encountered previously: the tri¬ 
angulations of the 2-disk (a,b), the projective plane (c), and the annulus (d). The 
roots are the 1-dimensional subcomplexes indicated with thick lines. In the first 
three cases j32(2f) = 0, so a root of X is simply a spanning tree of the codimension- 
1 skeleton. The annulus has j3i = 1, and its root is a unicyclic graph rather than a 
spanning 1-tree. 

In the case that {F,R) is a maximal rooted spanning forest, the root R was called 
a “relatively acyclic complex” in [DKM15], because the inclusion map ^ A in¬ 
duces isomorphisms //, (/?; M) —> Hi {X ; M) for all i <d, a condition equivalent to the 
vanishing of the relative homology groups for all i < d. Lyons [Lyo09] 




Fig. 4 Rooted trees and forests. 


used the term “cobase” for the complement of a root, referring to the fact that a 
/-cobase is a row basis for the coboundary map 

Definition 2.16. Let X be a c/-dimensional cell complex. A directed rooted span¬ 
ning forest of A is a triple {F,R, where {F,R) is a rooted spanning forest of X, 
and \Rd-i —is a bijection such that p C ^(p) for all /. That is, 

pairs each codimension-1 face p ^ R with a facet that contains it. 

When d = I, this formalism regards an oriented edge m —> v not as an ordered 
pair of vertices, but as the pairing of vertex u with edge uv. It is a standard fact that 
every rooting of a tree graph induces a unique orientation of its edges so that each 
edge is directed towards the root. In higher dimensions, a root can admit more than 
one orientation; this subject is discussed more fully in §5.3. 

Example 2.17. The rooted trees (a), (b) and (d) in Figure 4 each induce a unique 
orientation. The rooting of shown in (c) induces two distinct orientations: for 
example, edge 23 can be paired with either triangle 236 or triangle 234. See Figure 5. 
In general, the number of valid orientations is at least where Y is 

the relative complex {F,R), but equality need not hold, as shown by Bernard! and 
Klivans [BK]. 




Fig. 5 Two orientations of the rooted spanning tree of shown in Fig. 4 (c). 
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2.5 Other formulations of higher-dimensional trees 

The definition of a c/-forest fits naturally into matroid theory (for a general reference, 
see, e.g., [0x111]). The cellular matroid of a cell complex X is the matroid 
Mci{X) represented over M by the columns of dd{X). This is an extension of the 
definition of a simplicial matroid', see [CL87]. Proposition 2.13 then asserts that 
the spanning rZ-forests of X are precisely the independent sets in Md{X), and the 
maximal spanning forests are its bases. For a graph, this approach would define 
spanning trees as the bases of the graphical matroid. 

Higher-dimensional trees originate in the work of Bolker [Bol76], Kalai [Kal83] 
and Adin [Adi92], who studied spanning trees of skeletons of simplices and joins of 
0-dimensional complexes, respectively. The foregoing definitions, in varying gener¬ 
ality and essentially equivalent but for terminological differences, were formulated 
by Duval, Klivans and Martin [DKM09, DKMll, DKM15], Lyons [Lyo09], and 
Petersson [Pet09]. The present definition of a maximal spanning forest of a cell 
complex X matches the definition of “spanning tree” given by Petersson [Pet09, 
Defn. 3.13]. In addition, Petersson [Pet09, Defn. 3.14] defined a pure connected 
complex to be a tree if its top homology vanishes (i.e., if it is a forest accord¬ 
ing to Definition 2.5); by contrast, here Definition 2.6 emphasizes that the higher¬ 
dimensional analogue of connectedness is the vanishing of codimension-1 homol¬ 
ogy. Lyons [Lyo09] did not use the terms “tree” and “forest,” instead using matroid 
terminology. 

The literature contains many other generalizations of trees to higher dimension: 
see, e.g., [BP71, Dew74, Far02, HLM99, MV02, BPT15]. In general these do not 
have the same topological properties as the cellular trees considered here, but are 
defined so as to extend other useful properties of trees. For instance, Faridi’s defi¬ 
nition [Far02] extends the recursive description of a tree as a graph formed from a 
smaller tree by attaching a leaf. Faridi’s trees are simplicial complexes whose dual 
graphs are trees (such as the first, but not the second, tree in Figure 1); the presence 
of a leaf is important in the study of their facet ideals. In another direction, Mas- 
baum and Vaintrob [MV02] studied 3-uniform hypergraphs (which can be regarded 
as pure 2-dimensional simplicial complexes) and found a Pfaffian matrix-tree theo¬ 
rem enumerating their spanning trees. 


3 Enumerating spanning trees and forests 

This section will be somewhat more technical and algebraic, with the goal of giving 
precise statements of the various forms of the cellular matrix-tree theorem that have 
appeared in the literature and explaining how they are interrelated. 
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3.1 The classical matrix-tree theorem 


Let G = {y,E) be a loopless graph with n = \V\ and m= iLj. Let e(v,w) denote 
the number of edges between vertices i and j, let degg(/) be the degree of i, and 
let d be the signed incidence matrix of G (choosing an arbitrary orientation). The 
Laplacian matrix of G is L(G) = dd^. Explicitly, L{G) = where 


fdegoW 


if i = j, 

if i ^ j. 


The matrix-tree theorem enumerates spanning trees of a graph using its Laplacian 
matrix. 


Theorem 3.1 (Matrix-Tree Theorem). Let G be a connected graph, let ^{G) be 
its set of spanning trees, and let t(G) = \fi'(G}\. Then 

t(G) =- (6) 

n 

where the Xi are the non-zero eigenvalues ofL{G). Alternatively, for any i, j, 

t(G) = detZ (7) 

where L is the reduced Laplacian obtained by deleting the row and column for 
some vertex i (the choice ofi does not matter). 

The theorem can be proved by expanding the Laplacian or reduced Laplacian us¬ 
ing the Binet-Cauchy formula and observing that the non-vanishing terms corre¬ 
spond exactly to spanning trees. It is also possible to prove (7) using the dele¬ 
tion/contraction recurrence for t(G). 

Corollary 3.2. Let G be any graph (not necessarily connected), let fX (G) be its set 
of maximal spanning forests, and let t(G) = \ .^X(G)\. Suppose that G has nonzero 
Laplacian eigenvalues Xi,..., Xn-c and connected components Gi,..., Ge of orders 
n\,... ,nc. Then 

t(G) = = detL5 (8) 

ni---nc 

where Ls is the (n — c) x (n — c) matrix obtained by choosing a vertex in each Gi 
and deleting the corresponding c rows and columns from L. Moreover, the number 
of rooted forests ofG is simply the product of the nonzero eigenvalues ofL{G). 

Indeed, L(G) is the block-diagonal sum L(Gi) 0 • • • 0L(Gc), so its spectrum is the 
multiset union of their spectra. Meanwhile, a maximal spanning forest of L(G) is 
just the union of a spanning tree of each component, and the denominator in (8) 
is the number of possible choices of roots (since a root consists of a vertex chosen 
from each component of the graph). 
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Let w = {wg} be a family of commuting indeterminates indexed by the edges of 
K„. Let G be a graph on [n] and set Wg = 0 if e ^ ^(G). The weighted Laplacian 
matrix of G is L(G;w) = where 




if i = j, 

if i ^ j. 


The matrix-tree theorem may be extended to the weighted case, producing poly¬ 
nomial tree enumerators that carry finer combinatorial meaning. For example, let 
xi,...,Xn be indeterminates corresponding to the vertices of G and set Wij = xixj. 
Then the “Cayley-Priifer” weighted enumeration formula is 

Y. -(9) 

TeS’(Kn)i='^ 


3.2 Extending the matrix-tree theorem to arbitrary cell complexes 

The maximal spanning forests of a cell complex can be enumerated using a higher¬ 
dimensional analogue of the matrix-tree theorem. As described above, this idea orig¬ 
inated with Kalai [Kal83] and Adin [Adi92] and was extended further by others 
[CCK12, DKM09, DKMll, DKM15, Lyo09, Pet09]. 

Throughout, let X he a cZ-dimensional pure cell complex. Recall that the 
boundary map —)■ may be regarded as the incidence matrix from 

A:-cells to {k— l)-cells. Its adjoint map is the A:* coboundary map d^, which may 
be regarded as the transpose [5^] = [dkV ■ The A:* up-down, down-up, and total 
Laplacians are respectively 

Lt^iX)=dld,, Lr(X) +Lf, 

all of which may be viewed either as linear operators on or equivalently as 
square matrices of size |24|. They are all real, diagonalizable, positive semi-definite 
matrices. The spectrum of is the multiset union of the spectra of L'lf{X) 

and Lp{X) (because the up-down and down-up Laplacians annihilate each other). 
Moreover, the spectra of L'^‘^{X) and j {X) coincide up to the multiplicity of the 
zero eigenvalue. Therefore, the spectra of any of the families {L“'^(2f)}, {L^“(2f)}, 
determines the spectra of the other two. The results here are primarily 
phrased in terms of the up-down Laplacian, so we abbreviate by Lk{X), or 

simply Lj^ when the complex X is clear from the context. 

Higher-dimensional tree and forest enumeration requires keeping track of tor¬ 
sion. Accordingly, let Tor(A) denote the torsion summand of a finitely generated 
abelian group A, and define 


4(X) = |Tor(//i(X;Z))|. 
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We extend this notation: if S is any set of cells of X, we set tk{S) = \ Tor(/4(X5;Z))|, 
where Xg is the subcomplex generated by S. 

Let denote the set of maximal spanning A:-forests in X. As in (4) above, 

the fe* forest count of X is 


^k{X)= Y. tk{Tf. 

TeSkix) 

Note that when = 0 (i.e., the spanning forests are spanning trees), we 

have tk{T) = | Tor(f4(r;Z))| = \Hk{T-,'L) \ for each T. Moreover, when A = G is a 
graph, every summand is 1, so Ti (G) is just the number of maximal spanning 

forests. 


3.3 Weightings 

As for graphs, one can study weighted analogues of t(A) in order to obtain more 
refined information about its spanning trees and forests. 

Let w = {wf I F € X} be a family of commuting indeterminates indexed by the 
faces of X. The weighted forest enumerator of X is defined as 

t(A;w) = T,(A;w) = ^ tk-i{TfY\wF. (10) 

Tesr^(x) FeT 

When X is a simplicial complex, we frequently work with variables v, indexed by the 
vertices of X. and set Wa — IliGa for every face a. We will refer to this weighting 
as the vertex weighting for X, and denote it by v throughout. 

Let Dk denote the diagonal matrix with entries {wa : O’ € Xk). As in [CCK12, 
DKMl 1], one can define^ a “combinatorial” weighted boundary map dkO]/^, giving 
rise to the weighted Laplacian 


Lk-x{X-y/)^dkDkd*,. ( 11 ) 

Note that the combinatorial weighted boundary maps do not form a chain complex. 
Another approach [ADM, MMRW15, KL15b] is to define the weighted boundary 
map as j dkDj^ . These weighted boundaries do give rise to a chain complex, 
and to the “algebraic” weighted Laplacian: 

{X- w) = dkDk dl (12) 


* This weighted boundary was defined as d^Dk in [DKMll]. Here we adopt the equivalent con¬ 
vention of [CCK12] in order to eliminate squares on the indeterminates in weighted enumerators. 
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3.4 Cellular matrix-forest formulas 

The forest counts of a cell complex X can be computed from the spectra of its 
Laplacians. This relationship appears in several different forms in the literature, all 
of which share the same underlying principle. We refer to these results collectively 
as cellular matrix-forest theorems. 

The relationship between forest counts and the Laplacian becomes simplest when 
certain rational homology groups of X vanish. In [DKM09, DKMl 1 ] it was assumed 
that X was K-acydic in positive codimension (M-APC), i.e., that = 0 

for all k <d. This requirement is not necessary and so we weaken the hypotheses 
slightly. In fact, many complexes of combinatorial interest satisfy a stronger condi¬ 
tion: they are Z-acydic in positive codimension (Z-APC), i.e., Hk{X-,Z) =0 for all 
k < dimX. 

For S C Xji, the notation Ls indicates the square submatrix of Li^{X) restricted to 
rows and columns indexed by the cells in S. Define the pseudodeterminant X (M) 
of any square matrix M to be the product of its nonzero eigenvalues. 

Theorem 3.3 (Cellular-Matrix Forest Theorem). Let X‘^ be a cell complex with 
Hd-i{X-M)=Hd- 2 {X-M)=Q. LetR€ ^d-i{X) and S = Xd-i\Rd-i. Then 



(13a) 


(13b) 


Furthermore, for a general cell complex X‘^ with root R (see Definition 2.15), 

^d{X-,yv) = /%y'^L detT,(X;w). (13c) 

These three formulas are respectively [DKM09, Thm. 1.4], [DKMll, Thm. 2.8], 
and [DKM15, Prop. 3.5].^ In simple terms, (13a) and (13c) state that the forest 
counts of a complex can be computed from its reduced combinatorial Laplacian, 
and (13b) states that the forest counts can be computed via a product of eigenvalues 
of the combinatorial Laplacian. 

WhenX is Z-APC, the correction factors in Theorem 3.3 (e.g., td-i{X)) are equal 
to 1 and therefore may be omitted. 

Example 3.4. Let X be the equatorial bipyramid of Example 2.3; this complex is 
Z-APC. Let R be the root consisting of all edges that contain vertex 1. Then The¬ 
orem 3.3 states that the number of spanning trees equals the determinant of the 
5x5 reduced Laplacian formed by removing all rows and columns corresponding 
to edges in R-, this determinant is 15. 

^ There is a slight notational error in [DKM15], where the Lj- in Prop. 3.5 should refer to the 
reduced Laplacian formed by removing the rows and columns corresponding to cells of F. 



Simplicial and Cellular Trees 


19 


The formula (13b) can be iterated to obtain a closed-form expression for the tree 
count as an alternating product [DKMll, Cor. 2.10]. For simplicity, we state only 
the Z-APC case. 

Corollary 3.5. Suppose that = 0/or all k <d. Then 

( 14 ) 

i=0 


Example 3.6. Let X be the triangulation of the projective plane as given in Figure 3. 
The pseudodeterminants A, = A(L,) of its various Laplacians are: 

A 2 = 3‘^(3 + ^/5)^(3-^/5)^ = 3V, Al=6^ Ao = 6. 

Hence T2(2f) = A 2 A 0 /A 1 = 4. 


The right-hand side of Equation 14 is essentially the R-torsion (or Reidemeister 
torsion) of X with respect to the bases of the chain groups given by cells, in the 
sense of Milnor [Mil66, §8]. Catanzaro, Chernyak and Klein showed more gener¬ 
ally [CCK12, Thm. F] that the square of the R-torsion is the quantity 



covol Hk{X\'L) covol Bk{X\'L) 

w? 


^k+i{X) 


(-1/ 


where Bj^{X\'L) = (A). 

Catanzaro, Chernyak and Klein [CCK12] obtained a closely related determi¬ 
nant formula for Xci{X) in terms of the restriction Lb of Lii_i{X) to the space 
B = Bd-i = imdij-\{X) (which is invariant under the action of L^-i). For a lat¬ 
tice (discrete subgroup) .if C M”, define its covolume, covol.if, to be the volume 
of the parallelotope spanned by any (hence every) set of lattice generators. Then: 

Theorem 3.7 (Cellular Matrix-Forest Theorem, covolume expansion). [CCK12, 
Thms. C and D] 

Trf(X;w) = ^4z>i^detLB(X;w). 

(covol 

Theorems 3.3 and 3.7 are very close in spirit. Both express the torsion-weighted 
tree counts of X in terms of determinants of reduced Laplacians. Theorem 3.3 is 
useful in combinatorial applications, particularly where there is a natural choice of 
root R (such as when A is a shifted simplicial complex; see §4.3 below). Theorem 3.7 
is less combinatorial but has the advantage of expressing t(A) purely in terms of 
invariants of X, with no choice required. 

Lyons [Lyo09] gave another general formula for Tk{X) (which he denoted by 
hk-i{X)) without the requirements that Pd-iiX) = — 0. For a set of k- 

cells S C Xk, let tj^{S) denote the size of the group (ker^B^t)/((ker^Bjt nimQBjt+i) + 
ker2((9^|2f^y5)), and define a torsion-weighted enumerator for cobases by 
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hm= E 4(^^+i\5)4(5)^ (15) 

(A:+ 1)-cobases S 

Theorem 3.8 (Cellular Matrix-Forest Theorem, cohase expansion). [Lyo09, Prop. 6.2] 
Let S be a {d—\ )-cobase ofX, i.e., the complement of a {d — l)-root R. Then 

When = 0, these formulas reduce to (13b) and (13b) 

respectively. 

The algebraic weighting (12) can be used to obtain a slightly different matrix- 
forest formula. Versions of the following result appear in Aalipour, Duval and Mar¬ 
tin [ADM], Kook and Lee [KL15b, Thm. 4] and Martin, Maxwell, Reiner and Wil¬ 
son [MMRW15, Thm. 5.3]. 

Theorem 3.9 (Cellular Matrix-Forest Theorem with algebraic weighting). LetX 
be a cell complex of dimension d with weighting w, and suppose that i (A; M) = 

= 0- Let j (A) denote the algebraically weighted Laplacian ofX as 
defined in (12). Then 

The index d may be replaced with any k for which (A;IR) = = 0. 

Observe that setting all the weights w, to 1 recovers (13b). 

Just as iterating (13b) gives (14), iterating Theorem 3.9 yields an alternating 
product formula for weighted tree enumerators. For simplicity, we state only the 
Z-APC case. 

Corollary 3.10. Suppose that X is Z-APC of dimension d, i.e., HfX',^) = Q for all 
i < dimA. Then 

T,(A;w)= n 

(Tex ^=-1 

dim(7<i/ 

It is also possible to enumerate rooted forests in terms of products of Laplacian 
eigenvalues, as observed by Bernard! and Klivans. 

Theorem 3.11 (Rooted Cellular Matrix-Forest Theorem). [BK] Let X be a d- 

dimensional cell complex. Then the characteristic polynomial of the Laplacian 
Ld{X) gives a generating function for the rooted forests ofX: 



Simplicial and Cellular Trees 


21 


^ =Aci{La{X)+z-U). (18) 

(m 

Here Id is the identity matrix with rows and columns indexed by X^_ \. The result 
may also be extended to include indeterminate weights on the faces of the forest and 
root; see [BK, Thm. 18] and [MMRW15, Prop. 2.1], 


4 Enumeration of trees for specific complexes 

We next describe how the theorems of Section 3 have been used to obtain enumera- 
tive results for several families of cell complexes. In each case we have provided the 
full statements of the tree enumeration formulas but refer the reader to the references 
throughout for detailed treatments of these results. 


4.1 Skeletons of simplices 

Fix integers n> d>0, and let V = [n] = {1,2,..., n}. The complete cZ-dimensional 
complex on n vertices is the simplicial complex 

A = A„^d = {o’ C y I dimo < dj. 

Equivalently, A is the c/-dimensional skeleton of the simplex on vertices [nj. The 
complex A is Z-APC, with top Betti number 
Kalai [Kal83, Thm. 1] proved that 

(19) 

confirming Bolker’s [Bol76] guess that this formula was correct and generalizing 
Cayley’s formula t(K„) = n"^^. As described in §1.2, Kalai’s crucial observation 
was that the invariant is the right way to enumerate higher-dimensional trees. 
The actual number of spanning trees probably does not have a closed-form expres¬ 
sion; it deviates from Kalai’s formula first when n = 6 and d = 2, since the real 
projective plane can be triangulated with six vertices; see Example 2.12. To ob¬ 
tain formula (19), Kalai showed that the eigenvalues of the reduced Laplacian are 1 
(with multiplicity (^Ij)) and n (with multiplicity ("f) }, from which (19) follows 
immediately. 

The weighted Cayley-Priifer formula (9) also extends nicely to higher dimen¬ 
sions. Kalai [Kal83, Thm. 3'] proved that 

XdiAn^w) = +V2-\ -hv„)( d ) 


( 20 ) 
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where v is the vertex weighting. It is an open problem to give a bijective proof of 
this formula. It is not clear how to extend the Priifer code or Joyal’s bijection to 
arbitrary dimension, or how such a bijection should account for trees with torsion. 


4.2 Complete colorful complexes 

Let r and ni,...,be positive integers, and let Vi,..., Vr be pairwise-disjoint ver¬ 
tex sets with |y,j = Hi. The complete colorful complex is the (r— 1 )-dimensional 
complex defined as 

= {(Tcyiu---uy, |#((7ny,)< 1 v/}. 

The complexes are natural generalizations of complete bipartite graphs, 

which arise as the case r = 2. As described in §1.1, the history of higher-dimensional 
trees owes its origins to the appearance of complete colorful complexes in the trans¬ 
portation problems considered by Bolker [Bol76]. The pure full-dimensional sub¬ 
complexes of precisely the balanced simplicial complexes, which are 

significant in combinatorial commutative algebra; see, e.g., [Sta96, pp. 95-106]. 

As observed by Bolker, the complete colorful complex is Z-APC, with 

top Betti number (”i ~ !)■ Adin [Adi92, Thm. 1.5] gave the following closed 
formula for their tree counts. ForD C [r], set op = Y.ie[r\\D^i Kd = !)■ 

Then for all 0 < k < r — 1, 

/r-2-|0|\ 

(21) 

DC\r\ 

\D\<k 

In the special case k = r — 1, Adin’s formula becomes 

T.-i = n«; • (22) 

i=\ 

Bolker [Bol76, Thms. 27, 28] observed that (22) gives the exact count of trees if and 
only if at most two of the n,’s exceed 2 (because the smallest balanced torsion tree 
is a 9-vertex triangulation of the real projective plane 

Some additional special cases are worth mentioning. First, setting r = 2 in (22) 
recovers the formula (3) for the number of spanning trees of a complete bipartite 
graph. Second, if n, = 1 for all i, then is just the simplex on r vertices, so 

Tto = 0 for all D except D = 0, when nn) = 1, and a® = r, so (21) specializes to 
Kalai’s formula (19). Third, if n, = 2 for all i, then Kn^^...^nr is the boundary of the 
r-dimensional cross-polytope. Now = 1 for all D and op = 2(r— |D|), so (21) 
becomes 
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T*(^2,2,...,2) = n(2(r-£/))Q(''^V). 

d=0 

We will have more to say about this case in §4.6 below. 

Aalipour, Duval and Martin [ADM, Thm. 5.8] used the weighted matrix-tree 
theorem (Thm. 3.9) to prove a vertex-weighted version of Adin’s formula (21), as 
follows. Let V = {vjj I 7 € [r], f e [ny]} be a family of variables indexed by the 
vertices of For j G [r] and 0 < k < r — I, define 


ekj= E (-1)'"' 

r=l t=l DC\r\ teD 

\D\<k-\ 

Then 

/ r \ / \ 

TkiKn„...,n/^y) = Up?" n E *2 (23) 

V./=l / DClr] \je\r]\D J 

\D\<k 

Note that setting all Vj^t to 1 specializes pj to 1 and to Od, recovering (21). 

Complete colorful complexes play a significant role in the theory of cyclotomic 
matroids, as studied by Martin and Reiner [MR05]. In the irreducible case, is 
isomorphic to the dual of the simplicial matroid of the complete colorful complex 
K{n) := where n is the product of distinct primes pi,...,pr. That is, the 

bases of the matroid are precisely the complements of the facet sets of spanning trees 
of K{n). Musiker and Reiner [MR14] extended this work into a beautiful topological 
interpretation of the coefficients of the n‘^ cyclotomic polynomial In particular, 
Musiker and Reiner’s results explain why n = 105 = 3 • 5 • 7 is the first value for 
which the coefficients of are not all ± 1: it is the smallest n for which the complete 
colorful complex K{n) contains a spanning tree with torsion. In related work. Beck 
and Ho§ten [BH06] studied cyclotomic poly topes, which are dual to the generalized 
transportation polytopes considered by Bolker, and observed that patterns for their 
growth series break down for similar reasons involving torsion trees. 


4.3 Shifted simplicial complexes 

Let O’ = {to < • • • < is} and p = {70 < ” • < 7V} be sets of positive integers, listed in 
increasing order. Write o^p if oCp, or if s = r and /<, < ji^ for every k. This re¬ 
lation is a partial order, called the componentwise order. A simplicial complex A‘^ 
is called shifted if its vertices can be labeled so that in each dimension its 

faces form an order ideal with respect to this order: that is, if p G A and O ^ p, then 
O € A as well. Shifted complexes are an important class of extremal complexes; 
every simplicial complex can be transformed into a shifted complex in such a way 
that many combinatorial properties are preserved while the overall structure is sim- 
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plified. A signal success of the study of shifted complexes was the classification of 
all /-vector/Betti-number pairs for all simplicial complexes, by Bjorner and Kalai 
[BK88]. Shifted complexes of dimension 1 are known as threshold graphs and have 
generated much study in their own right; see [MP95]. 

The homology groups of a shifted complex A are free abelian; the Betti number 
j3jt(A) is the number of facets of A of dimension k not containing vertex 1. To see 
this, observe that contracting the subcomplex S = star 4 (l) produces a wedge of 
spheres, one for each facet not containing vertex 1. In particular, if A is pure then it 
is Z-APC, and 5 is a Z-acyclic spanning tree. For the remainder of this section, we 
assume that A is pure of dimension d on vertex set [n], 

Duval and Reiner [DR02] showed that shifted simplicial complexes are Lapla- 
cian integral. Specifically, let deg 4 (/) denote the number of facets containing ver¬ 
tex /, so that (deg^(l),...,deg 4 (n)) is a partition. Then the conjugate (transpose) 
partition is precisely the spectrum of the top-dimensional Laplacian (A), up to 
the multiplicity of the zero eigenvalue. This generalizes the corresponding result for 
threshold graphs, first observed by Hammer and Kelmans [HK96, Thm. 5.3] and 
stated in this form by Merris [Mer94]. 

Weighted enumerators for spanning trees of threshold graphs were calculated by 
Martin and Reiner [MR03, §6], and can be deduced as a special case of a more 
general bijective result of Remmel and Williamson [RW02, Thm. 2.4]. Weighted 
tree enumeration for pure shifted complexes was studied by Duval, Klivans and 
Martin [DKM09], as we now summarize. Introduce doubly-indexed indeterminates 
w = {^i,i I 0 < / < t/, 1 < y < n} and for each face a = {aq < • • • < ak], set 
Wa = nlo^i,ar In other words, Wa records not just which vertices occur in a, 
but their order within a. Theorem 1.5 of [DKM09] states that the eigenvalues of 
the weighted Laplacian Lii{A-,w) are polynomials in the weights, so that setting all 
weights to 1 gives an explicit proof of Laplacian integrality for shifted complexes. 
Specifically, the eigenvalues are in bijection with critical pairs of A: pairs (<y,p) 
of {k+ l)-sets of vertices such that a G A, p ^ A, and a is covered by p in the 
componentwise order. 

Let r = {p G A I l^p} and A={p€F|pU{l}€A} denote respectively the 
deletion and link of vertex 1 (which themselves are shifted complexes). Observe that 
for every critical pair (c7,p) = ((ao,.. .,«£/), {bo,.. ■,bd)), there is a unique index j 
for which bj = aj + l (and h, = a, for all i ^ j). Define the signature of (a, p) to be 
{S,T), where S = {aq, ..., Ay-i} and T = {1,2,..., A;}. Then [DKM09, Thm. 1.6] 


Trf(A;w) 



L wsuj' 
jgruji} 


(24) 


where Sig(F) denotes the set of signatures of critical pairs of F. The case d = I (i.e., 
A is a threshold graph) is equivalent to [MR03, Thm. 4]. Specializing Wij i—> vj 
for every i,j to obtain the vertex weighting v (called the “coarse weighting” in 
[DKM09]) produces the formula 
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I,(Zl;v) = f n E (») 

\j=2 ) V(s,T)GSig(r) iGru{i} / 


4.4 Color-shifted complexes 

As in §4.2, let K = be the complete colorful complex with vertex set Vi U 

• • • U Vr, where Vi = {v/j,..., A subcomplex A C /T is called color-shifted if 
it satisfies the following property: 

if C7 = {vpaj,V 2 ,a 2 ,---,Vr,aJ €2i and \<b<ai, then (7\{v,'fl,.}U{v,;i} e A. (26) 

Color-shifted complexes were introduced by Babson and Novik [BN06]. Struc¬ 
turally, color-shifted complexes are to complete colorful complexes as shifted com¬ 
plexes are to simplices. On the other hand, color-shifted complexes (even in dimen¬ 
sion 1) are not in general Laplacian integral. 

Let G be a color-shifted complex with r = 2, i.e., a graph. Assume that G is con¬ 
nected, and let n = m, m = n 2 . For !</<«, let Xi = max{y : vp,V 2 ,; G ^(G)}. 
Then X = (Ai,...,A„) is a partition, and it determines G completely: the edges 
of G correspond to the boxes in the Ferrers diagram of X, with Vi and V 2 corre¬ 
sponding to its rows and columns respectively. For this reason, 1-dimensional color- 
shifted complexes were studied under the name of Ferrers graphs by Ehrenborg and 
van Willigenburg, who found an elegant formula for their weighted tree enumera¬ 
tors [EvW04, Thm. 2.1]. Let G = G^^ be the Eerrers graph corresponding to a parti¬ 
tion X, and weight the vertices of Vi and V 2 by x = {xi,... ,x„} and y = {yi,... ,ym} 
respectively. Then 


T{G2,;x,y) =xi---x„yi---ym 0(3'! +•••+3'Ap-i) 11(^1 + '"+-^X„-i) 

p=l q=l 

where X is the conjugate partition to X. The proof of this formula in [EvW04] is 
recursive, proceeding by calculating T(G 2 ,;x,y)/T(G 2 ,';x,y), where X' is a partition 
obtained from X by removing a comer of the Eerrers diagram. Another proof, ex¬ 
pressing T(G 2 ,;x,y) in terms of the weighted tree enumerator of a related threshold 
graph, was given in [DKM09, §10.3]. 

Eor a color-shifted complex A with r colors, let v = be the vertex weight¬ 
ing, where q G [r] indexes the colors and i indexes the vertices of each color. 
Let A' be the subcomplex of A whose faces are {p € A | v,^.i ^ p for any q}. 
Note that every (r — 2)-face (7 G A is missing a unique color m = m{a)', set 
kia) = maxjy | aU {xm.t} € A}. Aalipour and Duval conjecture [AD] that 

T,_i(A;v) = PI {Xm,l+---+X„^kia)) 

q,i 


(28) 
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where e{q,i) =#{g G Ai—i \ Xqj G <7 and Xq/ j G <7 for some q' ^ q}. Aalipour and 
Duval have proven this formula in the case r = 3 (where a color-shifted complex 
corresponds to a plane partition). It reduces to (27) when r = 2, and to the weighted 
tree enumerator formula (23) when A is a complete colorful complex. 

Equation (28) resembles formulas (24) and (25) for weighted tree enumerators 
of shifted complexes. Specifically, the first factor on the right-hand side is a product 
of monomials indexed by faces with vertices numbered 1, and the second factor is a 
product of sums indexed by faces without vertices numbered 1. 


4.5 Matroid complexes 

Let M be a matroid on ground set V. The family J^{M) of independent sets of M 
forms a simplicial complex called a matroid complex. These complexes can be 
characterized in many ways; perhaps the simplest is that a simplicial complex is 
a matroid complex if and only if every induced subcomplex is pure. Matroid com¬ 
plexes are Cohen-Macaulay and Z-APC. (See Oxley [0x111] for a general reference 
on matroids, and Stanley [Sta96, Chap. Ill] for facts about matroid complexes.) 

Every matroid complex is Laplacian integral. This fact was proved first by Kook, 
Reiner and Stanton [KRSOO], and subsequently generalized by Denham [DenOl], 
who found an explicit eigenbasis for the Laplacian. An expression for t(j^(M)) 
can be extracted from Denham’s result (with difficulty) by applying (13b) or (14). 
One application of the Kook-Reiner-Stanton formula, given in [DKMll, §6.1], is 
an alternate proof of Adin’s formula (21), because the complete colorful complex 
Knj is the independence complex of a direct sum of rank-one uniform matroids. 

Kook and Lee [KL15b, KL15a] gave a simpler expression for weighted tree enu¬ 
merators of matroid complexes. Let a{M) — Tm{0,1), where Tm is the Tutte poly¬ 
nomial of M (see [B092]) and j3 (M) = (—1)''W ^. 4 C£(— (the Crapo beta 
invariant [Cra67], where r is the rank function; see also [Sta07, pp. 447^48]). Then 

/ X a(F)l3(M/F) 

eeE flats F \eeE\F / 

[KL15b, Thm. 9], where v is the vertex weighting. Consequently, 

T(j?f(M)) = ]^ (29) 

flats F 


4.6 Hypercubes 

The hypercube Qn is the space [0,1]” C M”, naturally regarded as a cell complex 
with cells {{0},{1}, [0,1]}". Note that the hypercube is not a simplicial complex. 
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but it has a tight combinatorial structure that enables enumeration of its spanning 
trees. 

The hypercube is contractible, hence Z-APC. The eigenvalues of L(){Qn,r) are 
0,2,... ,2j,...,2n, with multiplicities ("); this can be shown by directly construct¬ 
ing an eigenbasis as in [Sta99, pp. 61-62]. Applying formula (6) of the classical 
matrix-tree theorem yields the formula 

Ti(e„)=2-" Y[ 2|A|=22"-”-'flfe®. (30) 

AC[n] k=2 

This formula was discovered (in greater generality) by Cvetkovic [Cve71, (4.10)], 
using spectral techniques; see also [CDS80, pp. 75, 218]. The first purely combi¬ 
natorial proof was recently found by Bernard! [Berl2]. A generalization to higher 
dimension was obtained by Duval, Klivans and Martin in [DKMl 1, Thm. 3.4], who 
showed that the nonzero eigenvalues of Li^{Qn) are 2{k -f 1),.. ■,2{i -f 1),... ,2n, 
with multiplicities (^^^')(")- Applying the alternating-product formula (14) gives 
the tree count 

Tk{Qn) = n (31) 

i=k+i 

Weighted tree enumeration of hypercubes proceeds similarly. Introduce indeter- 
minates {qi,Xi,yi \ [n]} and weight each face/= (/i,...,/„) e {{0}, {1}, [0,1]}" 

by the monomial 


Wf = 



Now consider the algebraically weighted Laplacian w) (see (12)) for k>0. 

By [DKMll, Thm. 4.2], its nonzero eigenvalues are the rational functions u{A) = 
+ Qilyi^ indexed by subsets A C [n] with |A| > k, with multiplicities 
Applying the weighted alternating product formula^ (Corollary 3.10) yields 

T:k{Qn;y/) = {qvqn)^'‘=^-'^ - (u{A)Ylxjyj 

AC[n] \ jeA 

\A\>k 



as conjectured in [DKMll, Conj. 4.3]. For a complete proof of this formula, 
see [ADM]. The case k = \ gives the weighted enumerator for spanning trees of 
cube graphs [MR03, Thm. 3]. 


^ One needs the convention that w% = 1, so that Z.^f(2n) is just the 1x1 matrix with entry 

I.v&(Q„)Wv = Wj=l{Xj+yj). 
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4.7 Duality and self-duality 

Let be cell complexes of dimension d. We say that X and Y are dual if there 
is an inclusion-reversing bijection (7 i-A (7* between their face posets that induces 
a duality between their cellular chain complexes. On the level of chain complexes, 
this means that dx.k = i^rd-k+i every k, so^ = L‘f_j^{Y). Indeed, T CX 

is a spanning fe-forest in X if and only if = {(7* | (7 € X \ L} is a spanning k- 
forest in Y [DKMl 1, Prop. 6.1]. Duality carries over to the weighted setting: given 
a weighting w on X, define the dual weighting w* on F by w*((7*) = Then 

L'jf{X',w) = L^“(T;w*), and likewise for the algebraically weighted Laplacians 
(see §3.3). By the weighted cellular matrix-tree theorem (Thm. 3.3), TkiX'.'w) = 
Td-k{Y for all k. As an application, the (n — l)-skeleton of the hypercube Qn is 
dual to the complete colorful complex /f 2 , 2,...,2 (with n 2’s), so the formulas (31) 
and (32) for spanning tree enumerators of Qn are dual to special cases of (21) 
and (23) respectively. 

A central reflex is an embedding of a self-dual planar graph G on the sphere 
so that the isomorphism G ^ G* (where G denotes the planar dual) is realized by 
the antipodal map on S^. Tutte [Tut79] proved that if G is a central reflex, then t(G) 
is the square of the number of self-dual spanning trees of G. Maxwell [Max09] 
generalized Tutte’s result to higher dimension. Say that a matroid on ground set 
E = {ei ,..., ei ,..., e„} is involutively self-dual if it can be represented by an 
n X 2n integer matrix M of the form [/|A^], where I is the nxn identity matrix and N 
is skew-symmetric (so that the fixed-point free involution E —>■ E given by e, ^e, 
gives an isomorphism M —>M*). If A is the independence complex of such a matroid 
(so dim A = d = n — 1), then [Max09, Thm. 1.3] 

tM)= E td-iiTf=( E • (33) 

TeSj{A] \self-duairG5rf{4) / 

As an application, let k be a positive integer and let A be a simplex on vertex 
set [2k-1-2]. Then Tk{A 2 k+ 2 ,k) = (2k-|-2)(*:) by Kalai’s formula (19). The number 
is even for all k > 0, so this tree enumerator is a perfect square. The simplicial 
matroid of the k-skeleton of A is involutively self-dual via the map (7 i--> [2k -f 2] \ (7 
so this perfect square phenomenon is explained by (33). 

On the level of weighted tree enumerators, the formula for Tk{^2k+2,k'-’'^) given 
by (20) is a perfect square when n = 2k + 2 and k is odd (since in that case (^" j) 
must be even). Maxwell [Max09, Thm. 1.8] proved that this perfect square phe¬ 
nomenon is combinatorial, in the sense that the square root of Xk{A 2 k+ 2 ,k\'^^) has 
the combinatorial interpretation one would expect from (33): 


In [DKMll], it was erroneously stated that L'f{X) = In fact, duality interchanges up- 

down and down-up Laplacians. 
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Y, = Y tk-i{T)wT 

TeSkiA) \self-dual 765^(4) 

The proof identifies the left-hand side with the determinant of a skew-symmetric 
matrix, and the weighted enumerator for self-dual trees with its Pfaffian. 

Self-dual balls of even dimension were considered by Martin, Maxwell, Reiner 
and Wilson in [MMRW15]. If X is a self-dual cellular ball of dimension d = 2k 
and d = di^{X), then A(L<;_i(X)) = More generally, con¬ 

sider a weighting w such that Wa = Wo-*, where a* is the dual cell to a. Then 
A(L^*®j(X;w)) = T/t(X; w)^, a consequence of [MMRW15, Cor. 6.10]. If in addition 
the self-duality on X is induced by an antipodal map, then it can be shown that X 
has an orientation giving rise to a symmetric or skew-symmetric boundary matrix d 
(depending on the parity of k). In either case = X{dd*) = X{d)^, so in fact 

= A((9) and more generally Ti^{X;w) = X{Dj^d) by [MMRW15, 

Thm. 7.4]. 


5 Connections, further directions and open problems 
5.1 Critical groups of cell complexes 

The theory of critical groups emerges from a discrete dynamical system that has 
been studied in mathematics, physics, probability, and even economics in several 
closely related forms, including the abelian sandpile model, the chip-firing game, 
and the dollar game. Original sources include [BTW88, Big99, BLS91, Dha90, 
Lor91]. See [LPIO] for a short accessible introduction to the subject, and [Mer05] 
for an overview. 

The critical group can be defined purely algebraically from the Laplacian of a 
graph. In higher dimension, the algebraic definitions have natural analogues using 
the tools developed earlier in the chapter. 

Let G = {V,E)he a graph and let L — L{G) be the graph Laplacian. The critical 
group K{G) is the quotient group 

K{G) = Tor(Z^/im L) = ker^o/imT — Z^'^l^'/im L 

where L is the reduced Laplacian (see Theorem 3.1). The isomorphism is due to 
Biggs [Big99, Thm. 4.2]. Regarding a vector c € as a configuration in which c, 
is the number of grains of sand (or chips, or dollars) on vertex v,, the elements of 
the critical group can be thought of as equivalence classes of configurations modulo 
a dispersive Laplacian action. A special system of coset representatives, the critical 
configurations, encodes the possible long-term stable behaviors of the system. 

By the matrix-tree theorem, the order of the critical group, and hence the number 
of distinct long-term configurations, both equal the number of spanning trees of G. 
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There are deep connections between critical configurations and spanning trees, as 
well as the combinatorics of parking functions, G-parking functions, and hyperplane 
arrangements; see, e.g., [BW97, CLB03, CP05, BCTIO, PPW13, BS13, GK15a]. 

Accordingly, the critical group Ki (A) of a cell complex X is defined as 

Ki{X) =Tox{I?^' / \mLi) = ke.xdi/imLi = /\mLs 

where Lj is the reduced Laplacian obtained by removing a torsion-free /-tree, 
provided one exists. The isomorphism was proven by Duval, Klivans and Mar¬ 
tin [DKM13, Thm. 3.4] in the simplicial case; the proof carries over to the more 
general cellular setting. In analogy to the graph case, the size of Ki{X) is the forest 
count Ti+i (X). 

These higher critical groups model discrete dynamical systems that generalize 
the sandpile model [DKM13, §5]. A state of the system described by Ki{X) con¬ 
sists of a number for each /-cell; this quantity can be regarded as “charge” (/ = 0), 
“current” (/ = 1), “circulation” (/ = 2), etc. The Laplacian action redistributes this 
quantity between /-cells that lie in a common (/ + l)-cell. When / = 0, this system 
is the sandpile model described briefly above. When / = 1, each state assigns a cur¬ 
rent to each edge, and redistribution takes the form of circulation around 2-cells. 
Because the configurations are contained in ker(9,, the /-dimensional flow is con¬ 
servative: the sum of oriented flow at every /-dimensional face containing a fixed 
(/— 1)-dimensional face is zero. When / = 1, the oriented flow on the edges neither 
accumulates nor depletes at any vertex. When / = 2, the face circulation at each edge 
adds up to zero. 

The model can be regarded as having a sink consisting of an /-dimensional span¬ 
ning tree. By specifying flow values on all /-faces that are not in the sink, the val¬ 
ues on the sink needed to make the flow conservative are uniquely defined; this is 
why the configurations can be described by the cokemel of the reduced Laplacian. 
Guzman and Klivans [GK15b] have recently identified meaningful systems of rep¬ 
resentatives of the higher critical group which reflect the long term behavior of the 
Laplacian action. 


5.2 Cut and Flow Spaces 

The critical group of a cell complex is closely related to its cuts and flows, which 
generalize the corresponding objects in graph theory. In a connected graph G, an 
(edge) cut is a collection of edges whose removal disconnects the graph, and a 
bond is a minimal cut. In matroid terms, a bond is a cocircuit: a minimal set of 
edges that meets every basis. The support of every row of the incidence matrix d is 
a cut, namely the set of all edges incident to a particular vertex, and more generally 
the cut space is defined as im(9*, that is, the row space of d. 

Meanwhile, a flow on G is an assignment of a “current” to each edge so that 
the total flow into each vertex equals the total flow out (called a conservative flow 
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above); that is, a flow is an element of the flow space ker(9 = (im(9*)^. In matroid 
terms, the minimal supports of flows are circuits of the graphic matroid (i.e., simple 
cycles in G). The cut lattice “^(G) and flow lattice -^{G) are defined in the same 
way as the cut and flow spaces, except regarding <9 as a map of Z-modules rather 
than vector spaces. 

Accordingly, for a cell complex the cut space, cut lattice, flow space, and 
flow lattice of X in dimension k are defined as 

Cuti(A) = imR d*, {X), 'rfiiX) = im^ d,*(X), 

Flow*. (X ) = kerR d;,(X), ^k(X) = ker^ d^(X). 

An k-dimensional bond of A is a minimal set of k-faces that support a non-zero 
vector in Cuti.(A), and a k-dimensional circuit of A is a minimal set of k-faces 
that support a non-zero vector of Flow<, (A). (These are exactly the cocircuits and 
circuits, respectively, of the cellular matroid of A; see §5.4 and [SWIO].) Flows and 
cuts can also be defined homologically: e.g., vectors in FloW(/(A) are just elements 
of Nd(X;R), while vectors in CuG(A) are sets of facets whose deletion increases 
codimension-1 homology. 

Given a spanning tree T of a connected graph G, the (signed) characteristic vec¬ 
tors of the fundamental bonds and fundamental cycles of T form bases of the cut 
and flow spaces respectively, and in fact they are Z-module bases for the cut and 
flow lattices. This is a standard fact of algebraic graph theory [GROl, Chap. 14] 
that carries over well to general cell complexes [DKM15, §§4—6]. Predictably, it is 
necessary to account for torsion in higher dimension. The entries of the characteris¬ 
tic vectors of the fundamental bonds and circuits of a cellular spanning tree T C A 
can be interpreted as cardinalities of certain torsion groups [DKM15, Thms. 4.11 
and 5.3], and in order to form Z-module bases of %(A) or ^^{X), the tree T must 
be Z-acyclic. 

For an integral lattice ^ (that is, a subgroup of Z"), the dual lattice is .if# = 
{w G M" I (w,v) G Z Vv G .if}, where (•,•) is the standard inner product. Thus 
.if C .if#, and .if#/.if is called the discriminant group (or determinantal group) 
of .if. Bacher, de la Harpe and Nagnibeda [BdlHN97] proved that the groups K{G), 
^#7'^, and are mutually isomorphic. (See also [GROl, 

Chap. 14].) The relationship between the cut and flow lattices, their discriminant 

groups, and the critical group of A is more complicated in arbitrary dimension. The 
isomorphisms are replaced with short exact sequences 

0 -G A(A)^’^#/^ ^ ^ 0, 

where £ is the “error term” Tor(//^_i (A;Z)) [DKM15, Thms. 7.6 and 7.7], and 

A* (A) is the cocritical group, defined roughly by dualizing the definition of the 
critical group [DKM15, Defn. 7.4]. In particular, the groups A(A), A*(A), ‘^#/‘^, 
and Z^(® are all isomorphic if and only if i (A; Z) is torsion- 

free. 
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Bacher, de la Harpe and Nagnibeda’s work in [BdlHN97] arose from viewing a 
graph as a combinatorial analogue of a Riemann surface, where the critical group ap¬ 
pears as the discrete analogue of the Picard group of divisors, or the Jacobian group 
of holomorphic forms. This point of view also informs the graphical Riemann-Roch 
theorem of Baker and Norine [BN07] and subsequent work connected to tropical 
geometry [Luoll, HMY12, HKN13]. An open problem posed in [DKM13, §6] is 
to understand how the combinatorics of critical and cocritical groups on cell com¬ 
plexes parallels the geometric theory of group invariants such as Chow groups on 
algebraic varieties. 


5.3 Trees and simplicial complex decomposition 

Recall from Definitions 2.15 and 2.16 that a directed rooted spanning forest of 
a cell complex is a rooted spanning forest {F,R), together with an orientation 

that pairs each {d — l)-face p outside the root R with a facet (j) G F that con¬ 
tains p. As seen in Example 2.17, rooted forests of dimension > 2 can admit more 
than one orientation: as shown by Bernardi and Klivans [BK], the number of valid 
orientations is at least \Hi{-i{F,R)\. Equality need not hold: for instance, rooted tri¬ 
angulations of the (torsion-free) dunce cap can admit multiple orientations. In an 
appropriate linear-algebraic sense, these excess orientations cancel in pairs to the 
size of the torsion subgroup. It is an open problem to give a combinatorial explana¬ 
tion of this cancellation; see [BK]. 

Iterating the construction of a directed rooted spanning forest gives a discrete 
(although not gradient) vector field in the sense of discrete Morse theory [Eor98, 
Eor02]. To be precise, let {F,R, 0) be a directed rooted (7-forest of X‘^. Then R is a 
{d — 1)-dimensional subcomplex of X‘^. Choose a directed rooted forest of R which 
further pairs faces of X‘^, this time of dimensions {d — 2) and {d— 1). Continuing 
this process, the vector field is then given by the union of the pairings, and the 
unmatched cells are enumerated by the Betti numbers of X. Eor example, consider 
the directed spanning 2-tree of shown in (a) or (b) of Eigure 5. This is the 

1-dimensional root, and the arrows indicate the bijection 0 : R —> \R. In the 
next iteration, choose a single vertex as the 0 -dimensional root, and orient all of the 
edges of R away from the root to match them with nonroot vertices. Einally, match 
the empty face with the root vertex. 

These recursive orientations arise in theorems and conjectures about decomposi¬ 
tions of simplicial complexes, with applications to their /- and /z-vectors. In partic¬ 
ular, certain results about decompositions of simplicial complexes can be usefully 
restated in the language of simplicial forests and orientations, as we now describe. 

Suppose that A is a simplicial complex that is acyclic over at least one coefficient 
ring (it does not matter which), so that the /-polynomial f{A,x) = is di¬ 

visible by X — 1. Stanley explained this algebraic fact combinatorially by showing 
that the face poset of A can be decomposed as A' 0 FI, where A' is a subcomplex 
and there is a bijection 77 : A' —> such that every (7 € A' is a maximal proper 
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subface of rj{(j) [Sta93, Thm. 1.2]. Duval [Duv94, Thru. 1.1] showed more gen¬ 
erally that every simplicial complex A admits a decomposition BO A' 0 £2, where 
B contains faces for each k and A\£2 are as just described. Duval’s result 

(which in fact holds for all regular cell complexes) led to a new proof of Bjorner 
and Kalai’s characterization of /-vector/Betti-number pairs of arbitrary simplicial 
complexes [BK88], without using algebraic shifting. In Stanley’s and Duval’s con¬ 
structions, the set is a flag of forests: an order Alter in the face poset of A such 
that for every k, the set 12^^ = {a € | dime = k} is the set of facets of a spanning 

A:-forest of A. 

More generally, Stanley conjectured [Sta93, Conj. 2.4] that if A is A:-uply acyclic, 
i.e., if the link of every face of dimension < k —2 is acyclic, then its face poset 
can be decomposed into Boolean algebras of rank k. The k = 1 case is Stanley’s 
result cited above. The general ce remains open, although Duval and Zhang proved 
a weakening of Stanley’s conjecture in [DZOl]. Stanley introduced this conjecture in 
the same context as the partitionability conjecture, which stated that every Cohen- 
Macaulay simplicial complex A can be partitioned into Boolean intervals whose 
tops are its facets, with ranks given by the h-vector. In light of Remark 2.14, there 
are implications 

every link in X has a 2f is Cohen- every link in X has 

Z-acyclic spanning tree Macaulay a spanning tree 

because by Reisner’s criterion [Rei76], a simplicial complex is Cohen-Macaulay if 
the link of every face is Z-APC. 

The present authors envisioned simplicial spanning trees as an inroad to the par¬ 
titionability conjecture, culminating in its disproof in [DGKM15]. Meanwhile, in a 
different direction, Katthan [Katl4] had discovered that the depth conjecture (pro¬ 
posed by Stanley [Sta82]; see also [PSFTY09, Herl3]) could be reduced to the 
study of certain lattices associated to simplicial spanning trees. While the disproof 
of the partitionability conjecture implies that the depth conjecture is also false in 
general, Katthan’s methods reveal a deeper connection between trees and combina¬ 
torial commutative algebra: e.g., he showed that spanning trees of simplex skeletons 
(which he called stoss complexes) arise as the Scarf complexes of certain extremal 
lattices [Katl4, Thm. 4.3]. 


5.4 Trees and Matroids 

As mentioned in §2.5, the collection of all spanning forests of a cell complex X 
forms the collection of bases of a matroid, the cellular matroid M{X), which gener¬ 
alizes the simplicial matroids of [CL87] and coincides with the well-known graphic 
matroid in dimension 1. Matroid structures such as cocircuits are useful in studying 
cellular trees (as in §5.2), so it is natural to ask what other matroid structures and 
invariants mean in the context of trees. 
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Krushkal and Renardy [KR14] studied a polynomial invariant on cell complexes 
that generalizes the Tutte polynomial of the cellular matroid. (For a general ref¬ 
erence on the Tutte polynomial, see [B092].) Let X be a cell complex whose 
geometric realization is an orientable manifold of even dimension. Its Krushkal- 
Renardy invariant Px{x,y,a,b) is a generating function for spanning c/-dimensional 
subcomplexes Y C X, in which the exponents of x,y record combinatorial infor¬ 
mation about the cellular matroid and those of a, b record topological information 
about the embedding of F in X. Setting a = b = I recovers the Tutte polynomial 
of M{X), and P satisfies a duality relation Px{x,y,a,b) — Px*{y,x,b,a), where X* is 
the dual cell complex [KR14, Thm. 14], generalizing the usual Tutte duality relation 
TM{x,y) = TM*{y,x) on matroids. 

Evaluating the Tutte polynomial TM{x,y) of a matroid M at {x,y) = (1,1) gives 
the number of bases of M. If M is the graphic matroid of a graph G then rM(i,i) = 
t(G), but this equality does not hold for most cellular matroids because of torsion. 
On the other hand, torsion can be modeled combinatorially via the theory of arith¬ 
metic matroids, introduced by d’Adderio and Moci [DM12, DM13b, DM13a] and 
further developed by Fink and Moci [FM12]. An arithmetic matroid is a matroid 
equipped with a multiplicity function m that models the number of torsion elements 
in an abelian group. The arithmetic Tutte polynomial studied in [DM12, DM13b] 
records this multiplicity. For a cellular matroid M(A^), the multiplicity of a set A 
of facets is |Tor(//rf_i(AUA<^_i;Z))p, and evaluating the arithmetic Tutte poly¬ 
nomial at {x,y) = (1,1) gives the torsion-weighted tree count of X. Bajo, Burdick 
and Chmutov [BBC 14] similarly modified the Krushkal-Renardy polynomial to in¬ 
clude a torsion factor, again obtaining the torsion-weighted tree count as a suitable 
evaluation of their polynomial invariant. 

The matroid structure for trees suggests that graph algorithms could carry over 
well to the general setting of cell complexes. Of particular interest is the problem 
of sampling random trees, proposed by Russell Lyons [Lyo09] and Igor Pak. In the 
graph case, there exist efficient algorithms for tree sampling. The best known is Wil¬ 
son’s loop-erased random walk algorithm [Wil96], which generates a uniformly dis¬ 
tributed random spanning tree. Extending Wilson’s algorithm to general cell com¬ 
plexes is likely to be hard, as there is no known efficient algorithm to sample random 
bases of an arbitrary matroid. (It is possible for some classes, including equitable 
matroids; see [May06]; Gorodezky and Pak [GP14] generalized Wilson’s algorithm 
to a restricted class of hypergraphs that do not coincide with cellular forests.) One 
could look for an algorithm that samples cellular trees either uniformly, or accord¬ 
ing to a distribution that includes torsion information. Indeed, this problem informs 
Lyons’ motivation in [Lyo09], in which he defined determinantal probability mea¬ 
sures on cellular trees based on the Laplacian. 
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